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We present a derivation of the Keldysh action of a general multi-channel time- dependent scatterer 
in the context of the Landauer-Biittiker approach. The action is a convenient building block in the 
theory of quantum transport. This action is shown to take a compact form that only involves 
the scattering matrix and reservoir Green functions. We derive two special cases of the general 
result, one valid when reservoirs are characterized by well-defined filling factors, the other when the 
scatterer connects two reservoirs. We illustrate its use by considering Full Counting Statistics and 
the Fermi Edge Singularity. 

PACS numbers: 73.23.-b, 73.50.Td, 05.40.-a 
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I. INTRODUCTION 

The pioneering work of Landauer— £ and Biittiker^ 
lay the foundations for what is now known as the scat- 
tering approach to electron transport. The basic tenet 
is that a coherent conductor is characterized by its scat- 
tering matrix. More precisely the transmission matrix 
defines a set of transparencies for the various channels or 
modes in which the electrons propagate through the con- 
ductor. As a consequence, conductance is the sum over 
transmission probabilities. Subsequently, it was discov- 
ered that the same transmission probabilities fully deter- 
mine the current noise, also outside equilibrium, where 
the fluctuation-dissipation theorem does not hold 5 -. 

Indeed, as the theory of Full Counting Statistics 6 -! 7 -! 8 - 
later revealed, the complete probability distribution for 
outcomes of a current measurement is entirely character- 
ized by the transmission probabilities of the conductor. 
The fact that the scattering formalism gives such an ele- 
gant and complete description inspired some to revisit es- 
tablished results. Thus for instance interacting problems 
such as the Fermi Edge Singularity*^ was recast in the 
language of the scattering approac h 11 ' 12 : 13 ' 14 . The scat- 
tering approach has further been employed successfully in 
problems where a coherent conductor interacts with other 
elements, including, but not restricted to, measuring de- 
vices and an electromagnetic environmen t 15 ' 16 ' 17 ' 18 . It 
is also widely applied to study transport in mesoscopic 
superconductors^. 

Many of these more advanced applications are unified 
through a method developed by Feynman and Vernon 
for characterizing the effect of one quantum system on 
another when they are coupled^ 6 -. The work of Feynman 
and Vernon dealt with the effect of a bath of oscillators 
coupled to a quantum system. It introduced the concept 
of a time-contour describing propagation first forwards 
then backwards in time. By using the path-integral for- 
malism, it was possible to characterize the bath by an 
"influence functional" that did not depend on the sys- 
tem that the bath was coupled to. This functional was 



treated non-perturbatively. A related development was 
due to KeldysbJi. While being a perturbative diagram- 
matic technique, it allowed for the treatment of general 
systems and shared the idea of a forward and backward 
time-contour with Feynman and Vernon. Applications 
involving the scattering approach require both the no- 
tion of the non-perturbative influence functional and the 
generality of Keldysh's formalism. Until now, the com- 
bination of the Feynman- Vernon method with the scat- 
tering approach was done on an case-specific basis. Only 
those elements relevant to the particular application un- 
der consideration were developed. In this paper we unify 
previous developments by deriving general formulas for 
the influence functional, or equivalently the Keldysh ac- 
tion of a general scatterer connected to charge reservoirs. 

The Keldysh action of a general scatterer can be con- 
sidered as a building block. Its "interface" is the set of 
fields x± CO ■ Through this interface, the actions of many 
conductors can be combined into quantum circuits. As 
in the case of classical electronics, a simple set of rules, 
applied at the nodes of such a circuit, suffice to describe 
the behavior of the whole networ k 23 ' 24 . 

The influence functional Z[x±] and the Keldysh action 
-4[x±] = li-^txi] depend on two sets of time-dependent 
fields X+(0 an d X-(0 corresponding to forward and 
backward evolution in time with different Hamiltoni- 
ans. Recently, the situation was considered where x± 
are time-independent but the scatterer was allowed to 
fluctuate in time 2 ^. We consider the case where also 
the fields x± are time-dependent. Functional derivatives 
with respect to these fields generate cumulants of the dis- 
tribution of outcomes for the measurement of the degrees 
of freedom coupled to x±- Since these fields enter the 
Hamiltonian of the scatterer as a time-dependent poten- 
tial energy term, their effect is captured by the scattering 
matrix. Since the fields differ for forward and backward 
evolution, the scattering matrices for forward and back- 
ward evolution differ. 

Our main result is summarized by a formula for this 
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Keldysh action. 

A[s] = Tr In 



1+G A-G 
— hs— — 



— Tr In s_. (1) 



In this formula, G is the Keldysh Green function char- 
acterizing the reservoirs connected to the scatterer 25 . 
It is to be viewed as an operator with kernel 
G(a,o/;c;t,if)S CtC > where the Keldysh indices a, a' G 
{+,— } refer to time-contour ordering, c, d E Z refer to 
channel space, and t, t' G R are continuous time indices. 
The dependence on the fields x± is carried by the time- 
dependent scattering matrix s, that also has Keldysh 
structure, owing to forward and backward time-evolution 
with different Hamiltonians. Explicitly, it is to be viewed 
as an operator with kernel s(a; c, c'; t)S a , a 'S(t — t') where 
indices carry the same meaning as in the kernel of G. 
This formula is completely general. 

1. It holds for multi-terminal devices with more than 
two reservoirs. 

2. It holds for devices such as Hall bars where parti- 
cles in a single chiral channel enter and leave the 
conductor at different reservoirs. 

3. It holds when reservoirs cannot be characterized by 
stationary filling factors. Reservoirs may be super- 
conducting, or contain "counting fields" coupling 
them to a dynamical electromagnetic environment 
or a measuring device. 

When the reservoirs can indeed be characterized by 
filling factors f(e), the Keldysh structure can explicitly 
be traced out to yield 



A[s+,sJ =Trln -/) + «+/ 



Tr In §-. 



(2) 



In this expression operators retain channel structure and 
time structure. The time-dependent scattering matrices 
s± have kernels s±(c,c';t)S(t — t') that depend on the 
field X±(t)- I n "time" representation, / is the Fourier 
transform to time of the reservoir filling factors, and as 
such has a kernel /(c; t, t')5 c _ c r diagonal in channel space 
and depending on two times. This formula is of the same 
type as the Levitov-Lesovik formula for zero frequency 
Full Counting Statistics (FCS) 8 , but contains informa- 
tion about finite frequencies due to the arbitrary time- 
dependence of s±. 

Another formula may be derived from Eq. |T]), valid 
for two terminal devices. Each terminal may still be con- 
nected to the scatterer by an arbitrary number of chan- 
nels. We denote the two terminals left (L) and right (R). 
In this case the reservoir Green function has the form 



G 



G L 






Gr 



(3) 



channel space 



where Gur) have no further channel space structure. 
Matrix structure in Keldysh and time indices (indicated 



by a check sign) is now retained in the trace, but the 
channel structure is traced out. Thus is obtained 



Trln 



1 + T, 



{G l [ X ±},Gr[x±}}~2 



(4) 

In this expression, the field dependence x± is shifted to 
the Keldysh Green functions Gl and Gr of the left and 
right reservoirs. This formula makes it explicit that the 
conductor is completely characterized by its transmission 
eigenvalues T n . 

The plan of the text is as follows. After making the 
necessary definitions, we derive Eq. ([1]) from a model 
Hamiltonian. The derivation makes use of contour or- 
dered Green functions and the Keldysh technique. Sub- 
sequently, we derive the special cases of Eq. ([2]) and Eq. 
dH) . While the formulas ^ and ((U have appeared in the 
literature before, as far as we know, there has not yet 
appeared a formal derivation. 

We conclude by applying the formulas to several 
generic set-ups, and verify that results agree with the 
existing literature. Particularly, we explain in detail how 
the present work is connected to the theory of Full Count- 
ing Statistics and to the scattering theory of the Fermi 
Edge Singularity. 



II. DERIVATION 
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FIG. 1: We consider a general scatterer connected to reser- 
voirs. The top figure is a diagram of one possible physical 
realization of a scatterer. Channels carry electrons towards 
and away from a scattering region (shaded dark gray) where 
inter-channel scattering takes place. Reservoirs are charac- 



terized by Keldysh Green functions Gi n 



(out) ■ 



These Green 



functions also carry a channel index, in order to account for, 
among other things, voltage biasing. In setups such as the 
the Quantum Hall experiment where there is a Hall voltage, 
Gi n will differ from Gout, while in an ordinary QPC, the two 
will be identical. The bottom figure shows how the physical 
setup is represented in our model. Channels are unfolded so 
that all electrons enter at and leave at z + . 



We consider a general scatterer connecting a set of 
charge reservoirs. We allow the scatterer to be time- 
dependent. A sufficient theoretical description is pro- 
vided by set of transport channels interrupted by a poten- 
tial that causes inter-channel scattering. We consider the 
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regime where the scattering matrix is energy-independent 
in the transport energy window. Since transport is purely 
determined by the scattering matrix, all models that pro- 
duce the same scattering matrix give identical results. 
Regardless of actual microscopic detail, we may there- 
fore conveniently take the Hamiltonian of the scatterer 
to be 



H =VF ^2 J dz V>£j (z) {~l5 m ,nd z + U m>n (z)} 1pn(z) 



Tires + Ht, 



(5) 



where 7i res represents the reservoirs, and 7ix takes ac- 
count of tunneling between the conductor and the reser- 
voirs. The scattering region and the reservoirs are spa- 



tially separated. This means that the scattering poten- 
tial u mn (z) is non-zero only in a region z~ < z < z + 
while tunneling between the reservoirs and the conduc- 
tor only takes place outside this region. Note that in our 
model, scattering channels have been "unfolded" , so that 
in stead of working with a channel that confines particles 
in the interval (— oo, 0] and allowing for propagation both 
in the positive and negative directions, we equivalently 
work with channels in which particles propagate along 
(—00,00), but only in the positive direction. Hence, to 
make contact with most physical setups, we consider — z 
and z to refer to the same physical position in a channel, 
but opposite propagation directions. 

We consider the generating functional 



J 



Z = e JK = Tr 



T + exp \ -i / dt H + (t) } p T~ expii dtW (t) 



(6) 



in which Ti.^ is obtained from TL by replacing u mn {z) 
with arbitrary time-dependent functions In 
this expressions T + exp and T~ exp respectively refer to 
time-ordered (i.e. largest time to the left) and anti-time- 
ordered (i.e. largest time to the right) exponentials. In 
the language of Feynman end Vernon^ this is known as 
the influence functional. It gives a complete character- 
ization of the effect that the electrons in the scatterer 
have on any quantum system that interact with. Fur- 
thermore, the functional Z generates expectation values 
of time-ordered products of operators as follows. Let Q 
be an operator 



f Z+ 

Q=/ J f dzipj n (z)q mn (z)ip n (z). 

ran z 

Choose u* n (z,t) = u mn (z) + x±{t)q mn {z) ■ Then 

(t- fn<?(*j)j r+ (n Q(t * 

M , x JV 



(7) 



Z[X}\ X=0 (8) 



By merging the power of the Keldysh formalism of 
contour-ordered Green functions with that of the Lan- 
dauer scattering formalism for quantum transport, we 
obtain an expression for Z in terms of the Keldysh Green 
functions in the reservoirs and the time dependent scat- 
tering matrices associated with u =t (z, t). 

The argument will proceed in the following steps: 

1. Firstly we introduce the key object that enables a 
systematic analysis of Z, namely the single particle 



Green function g of the conductor, 
equations of motion that g obeys. 



We state the 



2. We define the Keldysh action A — InZ, and con- 
sider its variation 8 A. We discover that 5A can be 
expressed in terms of g. 

3. We therefore determine g inside the scattering re- 
gion in terms of the scattering matrix of the con- 
ductor and its value at the edges of the scattering 
region, where the reservoirs impose boundary con- 
ditions. 

4. This allows us to express the variation of the action 
in terms of the reservoir Green functions G; n ( out ) 
and the scattering matrix s of the conductor. 

5. The variation 8 A is then integrated to find the ac- 
tion A and the generating functional Z. 



A. Preliminaries: Definition of the Green function 

The first step is to move from the Schrodinger picture 
to the Heisenberg picture. To shorten notation we define 
two time-evolution operators: 



U±{t f ,U) =T+exp {~*/ ' dt> w± ^')} 



(9) 



Associated with every Schrodinger picture operator we 
define two Heisenberg operators, one corresponding to 
evolution with each of the two Hamiltonians Ti.^. 



Q±(t) =U±(t f ,uy QU±(t f ,U 



(10) 
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In order to have the tools of the Keldysh formalism at our disposal, we need to define four Green functions 



9m,n i z i t\Z ,t) 
9m,n i z i Z ,t ) 
9m,ni z i t\ Z ,t) 
9m,n( Z ' z ; t ) 



-e A Tr 



e^Tr 



= e^Tr 



e^Tr 



U+ (ti , t ) ip m + (z, t)p ipl_ {z', t') (U- (ti , to)) * 
U+ (h , t )ipl + (z', t') Po ^ m - (z, t) (U-(h,t )y 
U+itu^poT- ^j f n _(z',t')^(z,t)^ (U-(h,to)y 



(11) 



Here the symbol T + orders operators with larger time ar- 
guments to the left. If permutation is required to obtain 
the time-ordered form, the product is multiplied with 
(— 1)™ where n is the parity of the permutation. Sim- 
ilarly, T~ anti-time-orders with the same permutation 
parity convention. 

The Green functions can be grouped into a matrix in 
Keldysh space 

j., J i'\ ( 9nun( z i z ,t ) gm.n( z i z ; ^ ) | 

9m , n{ z,t,z,t)- y g7 + n{Z)t . z ^ tl) g—( Zt *,zr,i!))- 

Notation can be further shortened by incorporating 
channel-indices into the matrix structure of the Green 
function, thereby defining an object g(z, t; z' , t'). The el- 
ement of g that is located on row m and column n, is the 
2x2 matrix g m ,n- 

The Green function satisfies the equation of motion 

{idt + VFid z — vfu(z, t)} g(z, t; z , t') 

- J dt"E(z; t - t")g{z, t"; z't') = S(t - t')8(z - z')l. 

(13) 

The delta-functions on the right of Eq. (fT5|) encode the 
fact that due to time-ordering g++ and have a step- 
structure 

— 9{z - z')5(t -t'- Z —^)8 mn + f(z, t; z't') (14) 
vp vf 

where / is continuous in all its arguments. The self- 
energy 

E(.;r) = -*%^%--z)-i%^( 2 -, + ) (15) 
Zt c Lt c 

results from the reservoirs and determines how the scat- 
tering channels are filled. It is a matrix in Keldysh space. 
The time t c is the characteristic time correlations sur- 
vive in the region of the conductor that is connected 
to the reservoirs, before the reservoirs scramble them. 
Gi n (out) ( T ) is the reservoir Green functions where elec- 
trons enter (leave) the scattering region, summed over 
reservoir levels and normalized to be dimcnsionlcss. This 



form of the self-energy can be derived from the following 
model for the reservoirs: We imagine every point z in a 
channel m outside (z~,z + ) to exchange electrons with 
an independent Fermion bath with a constant density of 
states v. The terms 7i ros and 7ix are explicitly 

H ICS = f dEv I dzEa} n (E,z)a m {E 7 z) 

Ht = ^ Cm I dEv j dzipl l (z)a m (E,z) 

+a m (E,z)ip m (z), (16) 

where the tunneling amplitude c m characterizes the cou- 
pling between the reservoir and channel m. More gen- 
eral reservoir models need not be considered, since, as 
we shall see shortly, the effect of the reservoirs is con- 
tained entirely in a boundary conditions on the Green 
function g inside the scatterer. This boundary condition 
does not depend on microscopic detail, but only on the 
reservoir Green functions G m ( out ) . 

We do not need to know the explicit form of the reser- 
voir Green functions yet. Rather the argument below re- 
lies exclusively on the property of G; n ( out ) that it squares 
to unity2£: 

j dt"G(t-t") in{out) G(t"~t') in{out) =S(t-t')l. (17) 
A differential equation similar to Eq. (TT3|) holds for gK 



B. Varying the action A. 

We are now ready to attack the generating functional 
Z. For our purposes, it is most convenient to consider 
A = InZ. We will call this object the action. Our 
strategy is as follows: We will obtain an expression for 
the variation SA resulting from a variation u(z,t) — > 
u(z,t) + Su(z,t) of the scattering potentials. This ex- 
pression will be in terms of the reservoir filling factors / 
and the scattering matrices associated with u(z,t). We 
then integrate to find A. 
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We start by writing 

5A=-iv F e A ^2 ( dtfdz (^, m (2,t) (*) - k",.n( z .t) W 



(18) 



where 

<V4(*)V>»(*)> + (<) 



= Tr 



T + enp\-i dt' H + {t')\^ m {z)^ n {z)r + ey^l-i I dt 1 H + {t')\ p Q T- exp Ji/ dt' H~ (t') 



WMM*))- (*) 



Tr 



T+ exp 



in 



dt' W+(t') f po^ - exp J i / dt' W-(t') I ^O^OjOT- exp 



in 



r 



(19) 



C. Expressing 8 A in terms of the Green function g. 

In terms of the defined Green functions, the variation 
SA becomes 



5A 



IVp 



E 



in 



dt l dz 



% (z,t-0 + \z,t) 

+ Su n,rn(z, t)g~ n {z, t + + ;Z, t)J 

= iv F I dt I dzTr[Su(z,t)g(z,t + O k ;z 7 t)] 



(20) 

The object Su is constructed by combining the chan- 
nel and Keldysh indices of the variation of the poten- 
tial. The trace is over both Keldysh and channel in- 
dices. The symbol k refers to the regularization ex- 
plicitly indicated in the first line, i.e. the first time ar- 
gument of g ++ (z,t — + ;z,t) is evaluated an infinitesi- 
mal time + > before the second argument, while in 
g~~ (z, t + 0~; z, t), the first time argument is evaluated 
an infinitesimal time + after the second. This is done 
so that the time ordering (anti-time ordering) operations 
give the order of creation and annihilation operators re- 
quired in Eq. ([15]). 

It proves very inconvenient to deal with the k regular- 
ization of Eq. (f2T)]) . It is preferable to have the first time 
arguments of both g ++ and g evaluated an infinitesi- 
mal time + before the second. Taking into account the 
step-structure of g ++ we have 

g(z, t + k ; z't') = g{z, t-0 + ; z', t') 

1- V 



+ —s(t-t'-- — -)i 



Here T3 is the third Pauli matrix 



1 

-1 



(21) 



acting in 



g(z, t — + ; z', t') for points z and z' inside the scattering 
region where u is non-zero, to the value of g at z~ where 
electrons enter the scatterer. For z < z' and t < t' , the 
equations of motion give 



g(z,t-\ 0'z',t-\ ) 

vf vf 

= s(z,t)g(z-,t-0 + ;z'-,t')^(z',t'), 



where 
s(z, t) — Z exp 



dz"u(z" , t 



vf 



(22) 



(23) 



The symbol Z indicates that the exponent is ordered 
along the z-axis, with the largest co-ordinate in the inte- 
grand to the left. Note that the potential u at position z 
is evaluated at the time instant t + (z — z~)/vf that an 
electron entering the scattering region at time t reaches 
z. Often the time-dependence of the potential is slow 
on the time-scale (z + — z~)/vp representing the time a 
transported electron spends in the scattering region and 
u(z, t + ) can be replaced with u(z, t). This is how- 
ever not required for the analysis that follows to be valid. 
Substitution into Eq. (|24|) yields 



(24) 



SA =v F \ dt Tr [w(t)g(z~,t - + ; z~,t)] 



dt lim S(t - t')Ti 



with 

w(t) = 



w(t)l 



dzs\z, t)5u(z, t 
s\t)Ss(t). 



-)s(z,t) 



(25) 



Keldysh space. The equations of motion allow us to relate 



In this equation z + is located where electrons leave the 
scatterer. Importantly, here Tr still denotes a trace over 
channel and Keldysh indices. We will later on redefine 
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the symbol to include also a trace over the (continuous) 
time index, at which point the second term in Eq. (|24p 
will (perhaps deceptively) look less offensive, but not yet. 
In the last line of Eq. ([25]). s(t) = s(z + ,t) is the (time- 
dependent) scattering matrix. We sent the boundaries 
to and t\ over which we integrate in the definition of the 
action, to — oo and oo respectively, which will allow us to 
Fourier transform to frequency in a moment. The action 
remains well-defined as long as the potentials u + and u~ 
only differ for a finite time. 



D. Relating g inside the scattering region to g at 
reservoirs. Imposing boundary conditions implied 
by reservoirs. 



There is no inhomogeneous term on the right-hand side, 
because we restrict z to be less than z~ . We thus find 



g(z -0 + ,e;z ,e') 



e ieAz/v F exp 



Gin(e) 
2L 



Az 



g{z 



Az,e;z 



(28) 



Here the correlation length l c is the correlation time r c 
multiplied by the Fermi velocity vf- Using the fact that 
G(e)i n squares to unity, it is easy to verify that 



Our task is now to find g(z~,t — + ;z~,t). Because of 
the t — t' dependence of the self-energy, it is convenient 
to transform to Fourier space, where 



g(z,e;z ,e r ) 

Gin (out)( e ) 



dtdt' e lEt g(z,t;z~,t')e~ 
dt e let G(t) in ( out ). 



-ie t 



(26) 



In frequency domain, the property that G in ( out ) squares 
to unity is expressed as G in ( out )(e) 2 = 1. (Due to the 
standard conventions for Fourier transforms, the matrix 
elements of the identity operator in energy domain is 
2ir5(e — e').) The equation of motion for z < z~ reads 



-ie + vf9 z + 



G in (e) 
2t c 



g{z,e;z ,e') = 0. 



(27) 



exp 



G in (e) 
2l c 

1 + G in (e) 



Az 



■ exp 



Az 
21. 



1-G in (e) (Az 
2 6XP 2L 



(29) 



Since spacial correlations decay beyond z~ , g{z~ — 
Az, e;z~,e') does not blow up as we make Az larger. 
From this wc derive the condition 



[1 + G ln (e)] g{z- -0+,e;z-,e') = 0. 



(30) 



Transformed back to the time-domain this reads 

dt" [5{t - t") + G in (< - t")] g(z- - 0+, t"; z~, t') = 0. 

We can play the same game at z + where particles leave the scatterer. The equation of motion reads 



— ie + vf9 z + 9(z — z + ) 
This has the general solution 
g(z,e;z',e') 



+ N G ou t(£) 



2t c 



g(z, e; z+, e') = 2ttS(z - z')6(e - e'). 



(31) 



(32) 



exp < ie 



[{z - z+)9{z - z+) - {z 1 - z+)6{z' - z- 



Gout(s) 

2L 



g(z' - 0+, e'; z', e') + —9(z - z')S(e - e') 

V F 



(33) 



We will need to relate the Green function evaluated at z < z + to the Green function evaluated at z > z + , and so we 
explicitly show the inhomogeneous term. The same kind of argument employed at z~ then yields the condition 



[l-G out (e)] 



](z + -0+,e;z + ,e') + —5(e 
v F 



0. 



(34) 



where the inhomogeneous term in the equation of motion is responsible for the delta-function. In time-domain this 
reads 



dt" [S(t-t")-G out (t-t")] 



](z + -0 + ,t";z + ,t') + —6(t"-t') 

VF 



= 0. 



(35) 
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It remains for us to relate g{z + — + , t + z z ; z + , t' + z z ) to <?(z — + , t; z , £'). This is done with the help 



of Eq. (|22|) , from which follows 



g(z + - + ,i+ - —;z+,t' + — ) = s(t)g(z~ - 0+, t; z~, i')s f (0- 



v F 



v F 



(36) 



We substitute this into Eq. {55]) , multiply from the right with s(f') and from the left with s'(t). If we define 
G' out (t,t') = s^(t)G out (t — t')s(t') the resulting boundary condition is 



dt" [S(t-t")-G' ont (t-t")] 



g(z~ - 0+ t": Z -,t') + —Sit 1 ' - t') 

VF 



= 0. 



(37) 



E. Finding the variation of the action in terms of 
the reservoir Green functions and the scattering 
matrix. 



where the trace is over time, channel and, in the first 
term, Keldysh indices. The operator s' is related to the 
scattering matrix s through s' = As. 



At this point, it is convenient to incorporate time into 
the matrix-structure of the objects G; n , Gout an d <?• The 
resulting matrices will be written without overbars. Thus 
for instance s will denote a matrix diagonal in time- 
indices, whose entry (£, t') is 6(t — t')s(t). Similarly the 
(*,*') entry of G in (out) is G in(out) (i - t'). Also let g~ be 
the matrix whose (t,t') entry is g(z~ — + ,t;z _ ,i'). In 
this notation Gf n = G out 2 = I and Eq. (J3TJ and Eq. 
(1371) read 



(I + G in )g- = 
(I - G out ) (g- + l/v F ) = 0. 



(38) 



These two equations determine g uniquely as follows: 
From the first of the two equations we have 



= G' 



out (/ + G in )g~ 



= -(I-G' out )g- + (I + G' out G in )g- 



(39) 



In the first term we can make the substitution — (/ 
G' out )g~ = (7 — G' out )/vp which follows from Eq. 
Thus we find 



vf I 



^out^i 



-(I-G' m 



Gh 



G OVL t + Gi, 



(40) 



and the last line follows from the fact that Gf n = G out = 
I. We have taken special care here to allow for different 
reservoir Green functions at z~ where particles enter the 
conductor and z + where they leave the conductor. In 
order to proceed we must now absorb the difference be- 
tween the two Green functions in the scattering matrix. 
We define A through the equation 



G 



A _1 GinA 



(41) 



and drop subscripts on the Green functions by setting 
G = Gj„. Substituted back into Eq. (pM)) for the variation 
of the action yields 



6 A = Tr 



5s'(l - G) 



1 

Gs' + s'G 



-Tr[«_(S_)t] , (42) 



F. Integrating the variation to find the action A. 

We now have to integrate 6 A to find A. This is most 
conveniently done by working in a basis where G is diag- 
onal. Since G 2 = 1, every eigenvalue of G is ±1. There- 
fore, there is a basis in which 



G = 



I 

-I 



In this representation s' can be written as 



(43) 



(44) 



Here the two indices of the subscript has the following 
meaning: The first refers to a left eigenspace of G, the 
second to a right eigenspace. A subscript 1 denotes the 
subspace of eigenstates of G with eigenvalue 1. A sub- 
script 2 refers to the subspace of eigenstates of G with 
eigenvalue — 1. In this representation, 

*-<»a?T«i-( o^r^- (45) 



so that 



5 A = Tr [8s' 22 (s^ 1 )'] ~ Tr [«_ (S_)t] , (46) 



and thus 



A = Trlns', 



22 - Tr m S 

e A = (Det s_) _1 Det s' 2 . 



(47) 



In these equations, s_ is the scattering matrix associated 
with TL~ as defined previously. Its time structure is to 
be included in the operations of taking the trace and 
determinant. 

Note that in the representation where G is diagonal, it 
holds that 



l + G ,1-G 



'12 
'22 



(48) 
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Due to the upper- (block)-triangular structure it holds 



that Det s' 22 = Det [±±^ 
result 



„/ l-G 



] leading to our main 



A = Tr In 



1 + G .l-G 



Tr In S- 



(49) 



where it has to be noted that many matrices have the 
same determinant as the above. Some obvious examples 
include 



I 

s' 22 

I 

s 21 s 22 



= (l + G)/2 + (l-G)s'(l-G)/4 
= (l + G)/2 + (l-G)s'/2. (50) 



will also assume that electrons enter and leave a channel 
from the same reservoir, so that G; n = G ou t and hence 
s' = s. We recall as well as that the Keldysh structure 
of the scattering matrix is 



s+ 
s_ 



(52) 



Here s± have channel and time (or equivalently energy) 
indices. s± is diagonal in time-indices, with the entries 
on the time- diagonal the time-dependent scattering ma- 
trices corresponding to evolution with the Hamiltonians 

n ± . 

With this structure in Keldysh space, we find 



III. TRACING OUT THE KELDYSH 
STRUCTURE 

Up to this point the only property of G that we relied 
on was the fact that it squares to identity. Hence the re- 
sult (Eq. |4"9")) holds in a setting that is more general than 
that of a scatterer connected to reservoirs characterized 
by filling factors. (The reservoirs may for instance be 
superconducting). In the specific case of reservoirs char- 
acterized by filling factors it holds that 



1- 2f(e) 2f(e) 

2- 2f(e) -l + 2/(e) 



(51) 



e A =Det 



x Det 



l + (s+-l)/ -(S+-l)f 

I_(l-/) + / 

1 



(53) 



We can remove the Keldysh structure from the determi- 
nant with the aid of the general formula 



Det 



q = Det(AD- ACA^B) 

= Dct(DA- CA^BA). (54) 



Here /(e) is diagonal in channel indices, and / m (e) is the 
filling factor in the reservoir connected to channel m. We 



Noting that in our case the matrices B and A commute, 
so that CA~ 1 BA = CB, we have 



e A = Det 



Det 



(Mi -/) + /) + 

-/)+«+/! Det (SI 1 ). 



- (Mi - /) + / - 1) (*+ - Dct 



(55) 



r 



IV. AN EXAMPLE: FULL COUNTING 
STATISTICS OF TRANSPORTED CHARGE. 

A determinant formula of this type appears in the liter- 
ature of Full Counting Statistics^ of transported charge. 
This formula can be stated as follows: the generat- 
ing function for transported charge through a conductor 
characterized by a scattering matrix s is 



Z[ X ] =Det 1+(SLX-1)/ 



(56) 



where s x is a scattering matrix, modified to depend 
on the counting field \ that, in this case, is time- 
independent. (The precise definition may be found be- 
low.) 



As a consistency check of our results, we apply our 
analysis to re-derive this formula. We will consider the 
most general setup, where every scattering channel is 
connected to a distinct voltage-biased terminal. To ad- 
dress the situation where leads connect several channels 
to the same terminal, the voltages and "counting fields" 
associated with channels in the same lead, are set equal. 

Let us start by defining number operators 

7V fe = fdz [9(-z -z) + e(z-z )}ipl(z)Mz) (57) 



The operator A4 counts the number of particles in chan- 
nel k that are located outside the interval [—Zq, Zq\. The 
coordinate zq is chosen to lie between the reservoirs the 
scattering region. With the full counting statistics of 
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transported charge, we mean the generating functional 
(a la Levitov) 

2(xi, ■■■,XN,t) 

J E fe XfeAT fc (0)/2 p -i £ fc XfeM, (t)J J2 k Xfc A4 (0)/2 



(58) 

The function Z(x, t) generates all moments of the joint 
distribution function for (nx, n®, . . . , njy) charges to be 
transported into terminal (1,2,..., N) in the time inter- 
val (0, t). Formally t is sent to infinity, and this causes a 
singularity in all irreducable moments of the distribution 
of thransported charge. Dealing with this singularity is 
a subtle issue, which we will not concern ourselves with. 
The interested reader is referred to the literatu re 26 ! 27 . 

The dynamics of the the operators A4(i) are deter- 
mined by the Hamiltonian 

' H = v f^2 dz ^™( z ) {~id z 8 m . n + u TOi „(z)} ip n (z). 

m,n 

As already mentioned, we choose to count charges outside 
the scattering potential region, so that every A4 com- 
mutes with the potential energy term in the Hamilto- 
nian. Using the short-hand notation ^ m XmNm = X-Ni 
we manipulate the definition of Z(x,t) to find 

Z{ x ,t) = (jx-X/2jHt e -i x .N e -mt e -ix.X/2 

= (e m ^e- m -^). (60) 
In this equation, the Hamiltonian TL X is defined as 



Vf^Jdz ip^) {-id z 6 m ,n + "fcl( z )} ^n{z) 



(61) 



The transformed potential is Um] n {z) — u m ^ n (z) + 
5 mtn ^ L (S(z ~ z ) — S(z + zq)). One way to verify this is 
to note that formally Ti x is related to Tt through a gauge 
transformation where the gauge field in each channel is 
proportional to 6{z— zq)+0(~ z~ Zq). The delta-functions 
in u&J arise as a gradient of the gauge field that appear 
in the transformation of the kinetic term in Ti.. 

The calculation of the full counting statistics has now 
been cast into the form of the trace of a density matrix 
after forward and backward time evolution controlled by 
different scattering potentials. Our result, Eq. (|55|) . is 
therefore applicable, with 



s± 



Zexp 



z + 



dz u (±x) (z) 



e TO/2 Soe ±,x/2 = s±x . 



(62) 

In this equation, x is a diagonal matrix in channel space, 
with entries § m ,nXm- Substitution into Eq. gives 



Z(x) = Dct [l + (stj x - 1)/ 
in agreement with the existing literature 8 . 



(63) 



V. TRACING OUT THE CHANNEL 
STRUCTURE. 

A large class of experiments and devices in the field of 
quantum transport is based on two terminal setups. In 
such a setup the channel space of the scatterer is natu- 
rally partitioned into a left and right set, each connected 
to its own reservoir. We are generally interested in trans- 
port between left and right as opposed to internal dy- 
namics on the left- or right-hand sides. The scattering 
matrices have the general structure 



s± = X 



r t' 
t r' 



X~ 



X = 



X} 



(64) 



Here r(r') describes left (right) to left (right) reflec- 
tion, while t (t 1 ) describes left (right) to right (left) trans- 
mission (t is not to be confused with time). These ma- 
trices have no time or Keldysh structure but still have 
matrix structure in the space of left or right channel in- 
dices. The operators X^(t) and X^(t) have diagonal 
Keldysh structure (denoted by the superscript ±) and 
diagonal time structure (here indicated by r to avoid con- 
fusion with the transmission matrix t). They do not have 
internal channel structure and as a result the Keldysh ac- 
tion is insensitive to the internal dynamics on the left- or 
right-hand sides. Our shorthand for the Keldysh scat- 
tering matrix will be XsX~ l where we remember that s 
has no Keldysh structure. 

We now consider the square of the generating func- 
tional Z and employ the first expression we obtained for 
it (Eq. |4"9"|) which retains Keldysh structure in the deter- 
minant. 



Dct 



1 + G 



XsX- 



l-G 



Dets f (65) 



Here we exploited the fact that s_ acts on half of Keldysh 
space together with the fact that s + = s_, i.e. s has no 
Keldysh structure, to write exp2Tr lns_ = Dets. We 
now shift X to act on G and define 

G = X- X GX P = Q = (66) 



The operators P and Q are complementary projection 
operators i.e. P 2 = P, Q 2 = Q, PQ = QP = and 
P + Q = I. Because of this, it holds that Det(P + sQ) — 
Det(P + Qs). Thus we find 



Z l = e 2A = Det(PsT + sQ) 



(67) 



The left channels are all connected to a single reservoir 
while the right channels are all connected to a different 
reservoir. This means that the reservoir Green function 
has channel space structure 



G 



G, 



G, 



(68) 
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where Gl and G R have no further channel space struc- 
ture. At this point it is worth explicitly stating the 
structure of operators carefully. In general, an opera- 
tor carries Keldysh indices, indices corresponding to left 
and right, channel indices within the left or right sets of 
channels, and time indices. However P, Q and s are di- 
agonal or even structureless, i.e. proportional to identity 
in some of these indices. Let us denote Keldysh indices 
with k, k' £ {+,—}, left and right with a, a' € {L,R}, 
channel indices within the left or right sets with c, d G Z 
and time t, t 1 € R. Then P has the explicit form 

P(k,k';a,a';c,c';t,t') = P{k,k';a;t,t')S a . a ,S c , d . (69) 



The projection operator Q has the same structure. The 
scattering matrix s has the structure 



s(fc, fc'; a, a'; c, c'; t, t') — s{a, a'; c, c')6k,k'8(t — t'). (70) 



We now use the formula Det I ^ ^ J = 

Det(A)Det(P - C A^ 1 B) to eliminate left-right structure 
from the determinant. 



Z z 



P L r^ + Q L r P L t ] + Qr? 
P R t' ] + Q L t P R r^ + Q R r' 

= Dot (P L r f + Q L r) Det [p fl r' f + Q R r' - (p R t' f + Q L t) [p l r^ 1 + Q L r~ l ) (P L < f + Qr*') 

= Dct(P L r t + Q L r) 



x Det 



Pfl(r' t -P L t' t rt h^ + ir 1 ~Q L tr-h')Q R -P R {P L t' ] 1 1' + Q L t'^ r~ 1 t')Q R \ (71) 



Here it is important to recognize that the reflection and 
transmission matrices commute with the projection oper- 
ators Pl,r and Ql,r- Furthermore, notice that, in term 
b, the projection operator P R always appears on the left 
of any product involving other projectors, while Q R al- 
ways appears on the right. This means that in the basis 
where 



Pr = 



I 




Qr = 




I 



(72) 



term b is the determinant of an upper block-diagonal ma- 
trix. As such, it only depends on the diagonal blocks, so 
that the term P R (. . -)Qr may be omitted. Hence 



b = Dot 



P R (r n -P L i'V-V) 

+ (r' - Q L tr~H')Qi 



(73) 



Now we invoke the so-called polar decomposition of the 
scattering matrh*2£ 



r = uyT - TV 
t = iv'VTu' 



t' = iu\/Tv 



r' = V'VT^TV 



(74) 



where u, u' , v and v' are unitary matrices and T is a 
diagonal matrix with the transmission probabilities T n 
on the diagonal. We evaluate term a in the basis where 



Pl and Ql are diagonal to find 

' f VT 



a = Det 



= Det 



Tit* 



Tu' 



(iVT 



-T 



(75) 



Where I = Pl + Ql = Pr + Ql is the identity operator 
I(k,k';c,c';t,t') = Sk,k'S c ,c'S{t — t') in Keldysh, channel 
and time indices. For term b we find 



b = Det 



P R VT^T + Pl 



T 



T 



, Qi 

Combining the expressions for a and b we find 
Z 2 = e 2A = Det [1 - T(PrQ l + PlQr)] ■ 



(76) 



(77) 



Using the fact that Pl(r.) = (1 + ^i(ii))/ 2 and Ql(r) 
: I — Gl(r))/2 and taking the logarithm we finally obtain 
the remarkable result 



A 



T 



l + ^{{G L ,G R }-2) 



(78) 



This formula was used in [15j to study the effects on 
transport of electromagnetic interactions among elec- 
trons. In [l7T | the same formula was employed to study 
the output of a two-level measuring device coupled to the 
radiation emitted by a QPC. 
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VI. FERMI EDGE SINGULARITY 

In this section we show how our formulas apply to 
a phenomenon known as the Fermi Edge Singularity. 
The system under consideration is one of the most el- 
ementary examples of an interacting electron system. 
The initial analysi o 9 ' 10 relied on diagrammatic techniques 
rather than the scattering approach or the Keldysh tech- 
nique, and was confined to equilibrium situations. Sev- 
eral decades later the problem was revisited in the con- 
text of the scattering approac h 13 ' 14 . An intuitive deriva- 
tion of a determinant formula was given. Here we ap- 
ply our approach to confirm the validity of this previous 
work. We find exact agreement. This highlights the fact 
that the determinant formulation of the FES problem is 
also valid for multi-channel devices out of equilibrium, an 
issue not explicitly addressed in the existing literature. 

The original problem^ 1 ^ was formulated for conduc- 
tion electrons with a small effective mass and valence 
electrons with a large effective mass, bombarded by x- 
rays. The x-rays knock one electron out of the valence 
band leaving behind an essentially stationary hole. Until 
the hole is refilled, it interacts through the coulomb inter- 
action with the conduction electrons. The x-ray absorp- 
tion rate is studied. Abanin and Levitov reformulated 
the problem in the context of quantum transport where 
an electron tunnels into or out of a small quantum dot 
that is side-coupled to a set of transport channels. 





FIG. 2: A schematic picture of the system considered. It 
consists of a charge qubit coupled to a QPC. The shape of the 
QPC constriction, and hence its scattering matrix, depends 
on the state of the qubit. A gate voltage controls the qubit 
level splitting e. There is a small tunneling rate 7 between 
qubit states. 



We prefer to consider a slightly simpler setup that ex- 
hibits the same physics. The setup is illustrated in Fig. 
([2]). A Quantum Point Contact (QPC) interacts with 
a charge qubit. The shape of the QPC constriction de- 
pends on the state of the qubit. The Hamiltonian for the 
system is 

H = Hi |1) {l\ + (H 2 +e) |2) <2|+ 7 (|1) (2| + |2) (1|) (79) 

The operators Hi (H2) describe the QPC electrons when 
the qubit is in state |1), (|2)). They differ by a potential 
energy term, describing the the pinching off of the QPC 



constriction depending on the state of the qubit. We 
may take both Hamiltonians to be of the form (Eq. [5|) 
that we wrote down for a general scatterer. The energy 
e is the qubit level splitting, an experimentally tunable 
parameter. The QPC may or may not be driven by a 
voltage bias V . 

QPC electrons do not interact directly with each other 
but rather with the qubit. This interaction is the only 
qubit relaxation mechanism included in our model. We 
work in the limit 7 — > where the inelastic transition 
rates r 12,21 between qubit states are small compared to 
the energies eV and e. In this case, the qubit switching 
events can be regarded as independent and incoherent. 

Now consider the qubit transition rate r 2 i, from state 
|1) to 1 2) as a function of the qubit level splitting e. To 
lowest order in the tunneling amplitude 7 it is given by 

,0 

T 2 i = 2 7 2 Re / dre leT lim exp^(r) 
J-00 

cxp^(t) = tr [e^e-^-^Voe - ^ 1 * ] . (80) 

This is the usual Fermi Golden Rule. The time t over 
which we integrate can be interpreted as the time when 
the qubit switches from |1) to |2). The trace is over 
QPC states, and po is the initial QPC density matrix. 
We see that the expression for T21 contains an instance 
of the Keldysh action A that we have calculated. The 
correspondence requires us to set 

n + (t) = n 1 + (n 2 -ni)e(t-T)e(-t) 
n-(t) = Hi. (8i) 

In order to conform to the conventions of the existing 
literature, we write Z in the form where the Keldysh 
structure has been removed (Eq. [4J: 



A(t) = tr In J_(l — /) + 1+ (t)/ - tr In s 



(82) 



In this formula, s_ is the scattering matrix correspond- 
ing to H~ = Hi when the qubit is in state |1). It is 
proportional to identity in time- indices. The scattering 
matrix s+(t) corresponds to H + . It is still diagonal in 
time-indices but the diagonal elements s+(r)t are time- 
dependent. If we take the time it takes an electron to 
traverse the conductor to be much shorter than other 
time-scales such as the attempt rate of charge transfers, 
then 



s+(T) t = h + (s 2 - §i)0(t - T)fl(-t) 



(83) 



where §2 is the scattering matrix associated with H2 
when the qubit is in state |2). This expression first ap- 
peared in [13[ . In the language of the original diagram- 
matic treatment of the FES problerr&i^, it represents the 
total closed loop contribution. 

We may also write this closed loop contribution as 



Det 1 + (s\s 2 - l)n(r)/ 



(84) 
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where II is a diagonal operator in time-domain with a 
kernel that is a double step function 



n{T) t ,t> = o{-t)Q{t-T). 



(85) 



and the scattering matrices s\ and §2 no longer have 
time-structure. We may work in the channel space basis 
where s\§2 is diagonal. Its eigenvalues are e lXk . Suppose 
we are in zero-temperature equilibrium, then the filling 
factor / is the same in every channel. In the fourier 
transformed energy basis / is simply a step function: 



= S(e - e')9(-e) 



Thus one finds 



JjDet 



i)n(r)/ 



(86) 



(87) 



This determinant contains no channel structure any 
more. Operators only have one set of indices (time, or 
after Fourier transform, energy). II is a projection op- 
erator, diagonal in time-domain while / is a projection 
operator in energy domain. Such a determinant is known 
as a Fredholm determinant. 

The resulting transition rate is^ i 10 i 13 



r 21 (e) = 0(- 



e> \e\ \E C . 



where E c , is a cut-off energy of the order of the Fermi en- 
ergy measured from the bottom of the conduction band. 
The exponent a is known as the orthogonality exponent. 
It may be calculated by evaluating the Fredholm determi- 
nant analytically with Wiener-Hopf method. It is given 
in terms of the scattering matrices a a 12 ' 13 



1 

4^2 



tr In 



(4 



•S'2 



(89) 



with the trace being over channel indices. Inspired by the 
work of Abanin and Levito v 13 ! 14 we considered the case 
where the QPC is driven by a voltage bias. The results 
of our study may be found in [l8|. 



VII. CONCLUSION 

In this paper we have derived several expressions for 
the Keldysh action A for a general multi-terminal, time- 
dependent scatterer. This object is defined as the (loga- 
rithm of the) trace of the density matrix of the scatterer 
after evolution forwards and backwards in time with dif- 
ferent Hamiltonians: 



Tr 



T 



exp 



dtH + (t) 



to 



p T exp 



dtH~(t) 



(90) 



Our main result is a compact formula for the action 
in terms of reservoir Green functions and the scatter- 
ing matrix of the scatterer (Eq. [J). We have shown 
how to perform the trace over Keldysh indices explicitly 
when reservoirs are characterized by filling factors. Thus 
we obtained a formula (Eq. ^ belonging to the same 
class as the Levitov-Lesovik counting statistics formula. 
We have also explicitly performed the trace over chan- 
nel indices for a two terminal scatterer (Eq. 2]). In this 
case we demonstrated that the Keldysh action only de- 
pends on the scattering matrix through the eigenvalues 
of the transmission matrix. To illustrate the utility of 
the Keldysh action, and confirm the correctness of our 
results, we considered Full Counting statistics and the 
Fermi Edge singularity. We found that our results agree 
with the existing literature. 



1 R. Landauer, IBM J. Res. Dev. 1 223 (1957). 

2 M. Biittiker, Y. Imry, R. Landauer and S. Pinhas, Phys. 
Rev. B 31 6207 1985. 

3 M. Biittiker, Phys. Rev. Lett. 65 2901 (1990). 

4 M. Biittiker, Phys. Re. B 46 12485 (1992). 

5 Ya. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 

6 L. S. Levitov and G. B. Lesovik, JETP Lett. 58 230 (1993). 

7 L. S. Levitov, H. -W. Lee and G. B. Lesovik, J. Math. 
Phys. 37 10 (1996). 

8 I. Klich, Quantum Noise in Mesoscopic Systems, ed. Yu. 
V. Nazarov, Kluwer (2003). 

9 G. D. Mahan, Phys. Rev. 163 612 (1967). 

10 P. Nozieres and C. T. De Dominicic, Phys. Rev. 178 1097- 
1107 (1969). 

11 K. A. Matveev and A. I. Larkin, Phys. Rev. B 46 15337 
(1992). 

12 K. Yamada and K. Yosida, Prog. Th. Phys. 68 1504 (1982). 

13 D. A. Abanin and L. S. Levitov, Phys. Rev. Lett. 93 Art. 



No. 126802 (2004). 

14 D. A. Abanin and L. S. Levitov, Phys. Rev. Lett. 94 Art. 
No. 186803 (2005). 

15 M. Kindermann and Yu. V. Nazarov, Phys. Rev. Let. 91 
Art. No. 136802, (2003). 

16 M. Kindermann, Yu. V. Nazarov and C.W.J. Beenakker, 
Phys. Rev. B 69 Art. No. 035336, (2004) 

17 J. Tobiska, J. Danon, I. Snyman and Yu. V. Nazarov, Phys. 
Rev. Let. 96 Art. No. 096801 (2006). 

18 I. Snyman and Yu. V. Nazarov, Phys. Rev. Lett. 99 Art. 
No. 096802 (2007). 

19 C. W. J. Beenakker, Transport Phenomena in Mesoscopic 
Systems, ed H. Fukuyama and T. Ando, Springer (1992). 

20 R. P. Feynman and F. L. Vernon, Ann. Phys. (NY.) 24, 
118, (1963). 

21 L. V. Keldysh, Zh. Eksp. Teor. Phys. 47 1515 (19 64). 

22 A. G. Abanov and D. A. Ivanov. larXiv:0709.2898l (20071 

23 Yu. V. Nazarov, Superlattices Microstruct. 25 1221 (1999). 



13 



Yu. V. Nazarov, Phys. Rev. Lett. 73 1420 (1994). arXiv:0705.0099 (2007). 

J. Rammer and H. Smith, Rev. Mod. Phys. 58 323 (1986). 28 H. U. Baranger and P. A. Mello, Phys. Rev. Lett. 73 142 
D. A. Ivanov, H. -W. Lee and L. S. Levitov, Phys. Rev. B (1994). 
56 6839 (1997). 

J. E. Avron, S. Bachman, G. M. Graf and I. Klich, 



